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Manipulating energy levels while controlling the electron localization is an essential step for many 
applications of confined systems. In this paper we demonstrate how to achieve electron localization 
and induce energy level oscillation in one-dimensional quantum systems by externally controlling the 
deformation of the system. From a practical point of view, the one-dimensional potentials can be 
realized using layered structures. In the analysis, we considered three different examples. The first 
one is a graded quantum well between confining infinite walls where the deformation is modeled by 
varying slightly the graded well. The second systems is a symmetric multiple quantum well between 
infinite walls under the effect of biasing voltage. The third system is a layered 2D hybrid perovskites 
where pressure is used to induce deformation. The calculations are conducted both numerically and 
analytically using the perturbation theory. It is shown that the obtained oscillations are associated 
with level avoided crossings and that the deformation results in changing the spatial localization of 
the electrons. 


I. INTRODUCTION 

Over the past decade, manipulating electrons in 

S aantum-confined system has progressed remarkably [l|, 
and thus allows developing devices for many apMica- 
tions like quantum computing ii, spectroscopy [hi, Q , 
ultra-fast and sensitive optoelectronics B,!!! , optical tun- 
ability and negative heat capacity [ll[. In such 

devices, it is essential to have well controlled manipu¬ 
lation means. Different manipulation mechanisms are 
used currently such as superconductor circuits mui 
and spins in solids [i3| and molecules However, 

the controllability of confined quantum systems in gen¬ 
eral and in particular for selected energy levels is still a 
major challenge [TsI, [TqI that has stimulating the explo¬ 
ration of many new concepts [Ml. One of the largely 
investigated control techniques relies on the utilization of 
Stiickelberg oscillations. By passing an avoided crossing 
twice, the interference due to the dynamical phase be¬ 
tween transitions can be either constructive or destruc¬ 
tive. This is known as Landau-Zener-Stiickelberg (LZS) 
interference [illlllll. Usually, the levels’ oscillation 
occurs in the time domain. However, Nori and co-workers 
showed that this is achievable as well in the space domain 
using spatially inhomogeneous magnetic field [26|. Fur¬ 
thermore, similar oscillatory behaviors originated from 
different external influences and not intended for LZS 
interferometry, have been reported in other types of sys¬ 
tems. For example, using self assembled quantum dots 
[23 under a strain field, quantum dots in micro cavities 
under terahertz laser excitation [28| or coupled semicon¬ 
ductor nano rings with varying inter-ring distance [29|. 
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In this paper we present practical methods to manip¬ 
ulate electron localization and to cause energy levels’ os¬ 
cillation in quantum confined systems using controlled 
deformations that can, in principle, be applied at room 
temperature. These methods can be used to realized LZS 
interferometry where the oscillation is caused by various 
means. By deforming the quantum system, the system 
potential energy is slightly altered and hence the states 
and their energies are changed accordingly. The practical 
application of these ideas depends on the magnitude of 
the level oscillation and displacement of the localization 
centers, and how these are coupled to the externally con¬ 
trolled deformations. We found that the control of the 
state localization is easily achievable with deformations. 
The control of the energy levels requires careful engineer¬ 
ing addressing the binding strength of the electron to 
particular sites that can be obtained using, for example, 
quantum dots. The assumed controlled deformations are 
applied either by voltage biasing or pressure across 2D 
layered systems such as graded semiconductor quantum 
wells or stacks of 2D hybrid perovskites [ol, [33| . 

Thus, the problem becomes a piece-wise constant poten¬ 
tial and hence it has an analytical solution (3, where 
the eigenfunctions, which are exponential and trigono¬ 
metric functions with continuous logarithmic derivative, 
and the eigenenergies are obtained as solutions of tran¬ 
scendental algebraic equations. Also, the problem can 
be solved numerically using various methods like finite 
difference (FDM) [^. finite element (FEM) [^, , 

and spectral methods |4QL Idlj . In this work, we use a 
flexible FDM to calculate the eigen-pairs. 

The numerical calculations show that energy level os¬ 
cillations are achievable and manipulatable in the studied 
systems by controlled deformations. Furthermore, the 
obtained oscillations are associated with avoided cross¬ 
ings. Thus, it can be used to realize LZS interferometry 
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FIG. 1. (top left) The studied single-well structure, (bottom 
left and right panels) Numerical calculation of the first 20 
eigenenergies for the example system (L = 15, a = 2, = 3). 

The eigenenergies less than zero are shown in the bottom left 
panel while those greater than zero are shown in the right 
panel. 


using alternative means other than the time-varying field. 
It is shown also that beside the oscillations, the deforma¬ 
tion results in changing the spatial localization of the 
electrons. This should provide means to manipulate the 
electrons and hence can be used for other applications 
of quantum-confined systems. The next section presents 
three different scenarios of controlled deformations, fol¬ 
lowed by a final section with our general conclusions. 


II. CONTROLLED DEFORMATIONS 


In order to explore possible scenarios in which local¬ 
ization/delocalization transitions occur via controlled de¬ 
formations and how are they linked to energy level oscil¬ 
lations we studied three different confined systems. The 
first one is a graded quantum well between infinite walls 
where the deformation is modeled by varying spatially 
the graded potential. The second studied systems is a 
symmetric multiple quantum well between infinite walls 
under the effect of biasing voltage. The third system is a 
layered 2D hybrid perovskites where the pressure is used 
to induce deformation. Tuning the properties through 
pressure has been applied for purely inorganic systems 
1^ . 1^ . However, the softer character of organic lay¬ 
ers should lead to more pronounced effects in 2D hybrid 
materials For the calculations, we use a flexible 

finite FDM for the calculations as aforementioned 
and atomic units are assumed throughout the paper. 


A. Graded quantum well between infinite walls 


The system analyzed in this subsection is a square well 
between two infinite walls separated by a distance L. So, 


the potential (Figure-[T]) is simply 

{ 0 if {)<x<doTd^a<x<L 

—U if d < X < d^ a , 

CX 3 otherwise 

( 1 ) 

where d, a and U are all positive constants and d-\-a < 
L. The potential deformation is modeled by varying d, 
which is the distance between the left infinite wall and 
the edge of the well. Namely, the deformation represents 
to displace the position of the potential well with respect 
to the confining walls. As an example of the effects we 
show some of the results obtained with L = 15, a = 
2, and U = 3. In this case, we obtain two states with 
energy levels below zero (shown in lower left of Figure-[T]) 
that are bounded stated to the small well. The energies 
above zero as a function of d are shown in the right panel 
of Figure-dl up to the 20th level. It is clear that these 
states are oscillating with d and that the amplitude of 
the oscillations depends not in a monotonic way with the 
order of the level. 


In order to have an analytical perspective that could 
shed light into the different contributions that result on 
the observed oscillatory behavior, the non degenerate 
eigenenergies (above zero in the used model) are cal¬ 
culated using perturbation theory (PT). In this case, 
the unperturbed structure is assumed to be an infinite 
walls potential and the perturbation is due to the small 
finite well. So, the Hamiltonian can be rewritten as 
H = T Vo where 


Voir) 


0 if 0 < X < L 

oo otherwise 


and 


Vi(r) = 


—U if d < X < d a 

0 otherwise 


( 3 ) 


The well-known eigenvalues of Hq (T + Vq) are = 
^ (t) ’ where m = 1, 2 ,..., is the sum of peaks 
and valleys of the corresponding eigenfunctions, which 
are 




sin {kmx) 


( 4 ) 


Here, m is used in this subsection as the counter for the 
states with energies above zero and should not be con¬ 
fused with n that counts all the states. With the excep¬ 
tion of this subsection, n is used to describe the energy 
levels throughout this paper. The first order correction 
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FIG. 2. The densities for different states as indicated and selected values if d. As the small well is displaced from left to right 
there is a concurrent displacement of the density that can be concentrated on either side depending on the state level being 
observed. 


to the energies is given by 

pcI-\-q, 

= (m|Vi|TO) =-e / 'iljl,{x)dx 
Jd 

a 1 . ( a\ ( (2dH-a)^ 

= —U — + (7-sin miT— cos rmr --- 

L mir \ LJ \ ^ J 

( 5 ) 

It is clear from the above equation that the corrections 
are oscillating with d at a frequency of 2m7r/I/, which 
matches the oscillations obtained by the numerical cal¬ 
culations especially for higher levels. Eq.-|5] also implies 
that the oscillation amplitude is proportional to U and is 
sinusoidal with m and a, providing additional means to 
design and control the oscillations. In particular, the am¬ 
plitude dependency on m explains the behavior observed 
on Figured] showing stronger oscillations at intermediate 
values of m. 

The localization of the electron, as described by the 
density p{x) = |'0(x)p, is also affected by the defor¬ 
mation. In Figure [2] we show the resulted densities for 
the first four positive eigenenergies, n =3, 4, 5, and 6, 
for a selection of different d values. There is a clear 
shifting from one side of the system to the other on 
the most probable position of the particle as a function 
of the system state. For example, for the third state. 




FIG. 3. (color online) (left) The energies for the third and 
forth eignestates vs. d. (middle) The density for the third 
eigenstate for d = 6 and d — 7. (right) The density for the 
firth eigenstate for d = 6 and d — 7. 


the particle is localized in the right side of the well for 
d < {L — a)I2. The localization is shifted toward the left 
side for d > {L — a)/2. This shifting effect persists up 
to higher states, although the side chosen by the system 
depends, among other factor, on the odd/even character 
of the quantum number. 

So far we have shown that energy levels’ oscillations 
(for energies above zero) are directly related to the sys¬ 
tem’s deformation and therefore providing a way to ex¬ 
ternally control the particular values of the energy spec¬ 
trum of a system. As aforementioned in the introduction 
section, many similar oscillations, due to other causes, 
were reported and in many of them the oscillations are 
associated with crossing avoidance. The oscillations ob¬ 
served in this example are also associated with density 
swapping between states. In Figure O the energies and 
densities of the third and forth eigenstates are shown. In 
the left panel, it is clear that there is a crossing avoidance 
at d = 6.5. This is associated with a density swapping 
where for d = 6, the third state is localized in the left of 
the well while the fourth state is localized in the right. 
This swapping happens in the reverse order for d = 7. 


B. A symmetric multiple quantum well between 
infinite walls deformed by biasing 


In this subsection and the next we consider two ex¬ 
amples that could be practically achieved with currently 
known materials and technologies, allowing the experi¬ 
mental exploration of controlled energy level oscillations 
and the associated localization/delocalization of the par¬ 
ticles. In this case, we start with a symmetric double-well 
system confined between two infinite walls. Then, a bias¬ 
ing voltage is applied to deform the original system into 
a tilted one. The unbiased potential, represented in the 
sketch of Figure S] is mathematically described by the 
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FIG. 5. (left) The energies for the and 12^^ eignestates 
vs. Vh- (middle) The density for the eigenstate for Vh = 
—0.2 /iV and Vh = 0.2 /iV. (right) The density for the 12^^ 
eigenstate for Vb — —0.2 /xV and Vh — 0.2 /xV. 


C. A layered 2D hybrid perovskites deformed by 
applying pressure 


FIG. 4. The sketch shows the double-well symmetric struc¬ 
ture. The left figure displays the energy levels as a function 
of the biasing voltage Vh and the induced oscillations. The 
different colors represent the different energy levels. 


following potential: 

{ 0 if \x\ < d OT d ^ a <\x\ < L 

—U if d <\x\ < d^ a , 

oo otherwise 

(6) 

where d, a and U are positive constants and d-\-a < L. 
Equivalent systems can be made although with finite po¬ 
tential walls. For example, using semiconductor multi¬ 
ple quantum wells [43, Q ^ double quantum dot sys¬ 
tem [ 45 I, nil can be used to realize a device having the 
qualitative features of Eq. Having this examples in 
mind, the parameters used in this analysis are within the 
ranges of practical values. We chose L = 200, a = 40 
and U = 0.02 that represent 10 nm, 2 nm and 0.54 eV, 
respectively. By applying a biasing voltage 14, the po¬ 
tential is tilted by an additional term of —VbX providing 
the externally controlled deformation. In our example, 
the voltage is varied between -0.15 and 0.15 mV. The re¬ 
sulting energies are shown in Eigure |4] as a function of 
14 , where it can be seen the emergence of an oscillat¬ 
ing pattern having similar characteristics to the example 
of the previous section. Eor example, the amplitude of 
the oscillation is a function of the energy level, as it was 
found by the PT analysis presented above. 

Analogously to the case of Subsection A, in this ex¬ 
ample the energy oscillation has an associated crossing 
avoidance and localization transition. In Eigure [5] we 
show the energies and densities corresponding to the 11^^ 
and 12^^ eigenstates. In the left panel, it is clear that 
there is a crossing avoidance around I 4 = 0 mV. The 
concurrent density swapping happens upon a very small 
change in the biasing: for Vb = —0.01 mV, the 11^^ state 
is localized in the left of the well while the 12^^ state is 
localized in the right. The picture is completely reversed 
for Vb = 0.01 mV. 


The second possible practical deformation method that 
we consider relies on the application of an external pres¬ 
sure. To realize this system we propose a layered 2D hy¬ 
brid perovskites. These structures are composed of alter¬ 
nating layers of hybrid 2D framework of octahedra inor¬ 
ganic and organic cations providing the wells. Moreover, 
the system could be intercalated between softer poly¬ 
meric layers allowing large deformations and increasing 
the degree of external control by the application of mod¬ 
erate pressure as compared with purely inorganic crys¬ 
tals 1471-1^. These materials can be analyzed as super¬ 
lattices, as shown by Even et al. [9| and the references 
therein. 

The most studied 2D hybrid materials are 
[RNH 3 ] 2 (CH 3 NH 3 )^.,Pb^l 3 „+i, with R represent¬ 
ing an organic group. In these materials, 
forms the octahedra that include the (CH 3 NH 3 ) cations 
forming the wells. On the other hand, [RNH 3 J 2 consti¬ 
tutes the barriers and by changing the organic group 
(R), the equilibrium separation between the wells can 
be tailored [9[. The depth of the well is estimated to 
be Ui = 0.9 eV [H^. As for the well width, it depends 
on the number of the octahedra layers. Eor m =1, 2, 3, 
and 4, the thicknesses a are 0.65, 1.25, 1.85, and 2.5 nm 
respectively 0 . 

In this subsection, we study the effect of 
the pressure on the energy levels on the 
[RNH 3 ] 2 (CH 3 NH 3 )„. 3 Pb^l 3 „+l quantum wells with 
1.25 nm barrier. A finite number of wells will be used 
where the whole system is then placed between two 
insulators with a barrier height of U 2 = 2.7 eV, as 
described in Eigure [ 6 ] for the case of two wells [9[. In our 
approximation, the external pressure only affects the 
confining external layers by changing their thickness by 
a factor p. The eigenvalues show the avoidance crossing 
characteristics as in the previous cases, although now 
as a function of the deformation parameter p. This is 
explicitly displayed in Eigure (3 in which we display 
two states corpsonding to the case of m = 3. As the 
system is deformed, the 14^^ energy level has the particle 
shifting its preferential position from the central region 
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trons. This effect is quite general, as shown by the the 
three described examples, and therefore it can be real¬ 
ized in different ways. The external controlling param¬ 
eter could be an applied voltage or pressure, and other 
procedures are certainly possible. Obtaining a practical 
way to experimentally control the electron localization 
would enable the manipulation of charges in quantum 
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FIG. 6. (top) Model potential to represent 2D layer per- 
ovskites confined between insulating walls, in this case for two 
layers creating two confining wells. The deformation affects 
the external confining layers of width dp. (bottom) Numerical 
calculation of the first 40 eigenenergies for a two wells system 
with m =1, 2, 3, and 4 as a function of the deformation factor. 

between the wells for p = 0.6, to the outermost layers 
for p = 1.4. The energy level immediately above, the 
situation is exactly the opposite. 

III. CONCLUSION 

In this paper, we have analyzed the interplay between 
the detailed dimensions of a confining potential, energy 
level oscillations and particle’s localization. It was shown 
that the obtained oscillations are associated with level 
avoided crossings and that the deformation induces a 
change in the preferential spatial localization of the elec- 


FIG. 7. (top) Eigenvalues corresponding the 14^^ and 15^^ 
energy levels for the case of two wells and m = 3 as function 
of the deformation factor (the pressure p). The red vertical 
dashed lines correspond to the the selected deformations for 
which the densities are plotted for the two states as function of 
distance, (bottom) The densities of the two states as function 
of distance where the dashed green lines are corresponding to 
the edges of the insulating walls while the dashed red lines 
are corresponding to the small wells edges. 

confined systems and provide guidelines for the design 
of devices. The approach that we have used to describe 
the is based on simple ID problems, but the consistent 
presence of this effect on multiple scenarios suggests that 
this phenomena is general and could be realized in 3D 
systems at the nanometer scale. 
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